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' Abstract 

' Gross-Pitaevskii and nonlinear Hartree equations are equations of non- 

linear Schrodinger type, which play an important role in the theory of 
Bose-Einstein condensation, Recent results of Aschenbacher et al |3 demon- 
strate, for a class of 3— dimensional models, that for large boson number 
(squared norm), M, the ground state does not have the symmetry 
^ ' properties as the ground state at small N. We present a detailed global 

, study of the symmetry breaking bifurcation for a 1— dimensional model 

fSJ ' Gross-Pitaevskii equation, in which the external potential (boson trap) is 

' an attractive double-well, consisting of two attractive Dirac delta func- 

0^ , tions concentrated at distinct points. Using dynamical systems methods, 

' we present a geometric analysis of the symmetry breaking bifurcation of 

' an asymmetric ground state and the exchange of dynamical stability from 

the symmetric branch to the asymmetric branch at the bifurcation point. 



> ■ 1 Introduction 

■ The experimental realization of Bose-Einstein condensation has given great im- 

d ' petus to the study of the Gross-Pitaevskii and nonlinear Hartree equations. 

These are equations of nonlinear Schrodinger (NLS) type, with a potential hav- 
ing both hnear and nonlinear parts. Such equations arise in the study of large 
dilute systems of bosons; A^-body quantum mechanics of coupled bosons, where 
N tends to infinity and the coupling strength tends to zero, |H IllL 151 IT]. Of 
fundamental interest are the properties of the ground state of such systems. 
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Because of the nonlincarity, the ground state can exhibit non-trivial transitions 
in structure. 

Indeed, that such a transition takes place was established in the recent paper 
They consider the nonlinear Hartree energy functional 

Hg[^,^^l [ \V'^ix)\^ +vix)\^{x)\^ +g\-^ix)mv^\^\^){x) dx. (1.1) 
^ Jtr'' 

Here, v{x) is an external potential (the boson trap), and V{x) is the two-body 
Coulomb interaction potential V{x) — \x\~^ between bosons. The ground state 
of a system of A/" bosons is characterized by solution to the minimization problem 

inf{i?,[vI/,^']:||vI/||2^AA}. (1.2) 

For spatial dimensions d > 2 and the case in which 17 < is fixed, corresponding 
to attractive interatomic forces (negative scattering length) , it is shown that for 
a class of potentials v{x), e.g., the double well potential, that if A/" is sufficiently 
large, then symmetry breaking in the ground state occurs, in the sense that the 
ground state of H1.2|) does not have the same symmetries as v{x). The proof 
in 3 is based upon the intuition in the double-well case, for example, that 
for small amplitude states (M small) the nonlinear potential is negligible, and 
therefore the ground state should be a small distortion of the linear ground state 
(g = 0), and have the same symmetries as the linear ground state. However, 
by careful construction of trial functions, one can see that for M sufficiently 
large, JV > A/thresh, it is energetically preferable for the ground state to be 
concentrated in one well or the other, but not equally in both. 

The size of A/thresh and the dynamical stability properties of the various 
coexisting states (symmetric and asymmetric) is unaddressed. In this work we 
address these questions in the context of the simpler model, the one-dimensional 
Gross-Pitaevskii equation (GP/NLS) 

i^t = ~ 2\i;\^^ + e vix)^ (1.3) 

where v{x) denotes the attractive double well potential 

v{x) ^ -S{x-L)-5{x + L), (1.4) 

with conserved Hamiltonian energy functional 

+v{xmx)\^-^mxr^ dx (1.5) 

and conserved particle (boson) number 



A/'[V',V5] = 




Symmetry breaking bifurcations have been observed in this model jBj and in 
the analogous double square well model [T2\ , by piecewise analytical solution in 



diplx) 



dx 
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terms of hyperbolic secant and Jacobi elliptic function solutions of the trans- 
lation invariant NLS equation (v = 0), followed by numerical solution of the 
nonlinear algebraic equations resulting from the required continuity and jump 
conditions at singular points of the equation's coefficients. 

We present a simple geometric analysis demonstrating 

(A) the symmetry breaking of the ground state, as a bifurcation of an asym- 
metric branch of solutions from the branch of symmetric solutions 

and 

(B) the exchange of dynamical stability from the symmetric branch to the 
asymmetric branch. In particular, beyond the bifurcation point the sym- 
metric state is linearly exponentially and nonlinearly unstable, whereas the 
asymmetric state is nonlinearly orbitally stable. 

We note that our analysis docs not require the specific cubic nonlinearity and is 
applicable for general local nonlincarities, where IV'PV' is replaced by /(IV'P)V'j 
with /(•) of a general class. 

A bifurcation diagram showing the symmetric, antisymmetric and asymmet- 
ric states is presented in figure ^ The solid curve corresponds to the nonlinearly 
stable ground state. As A is increased (— is the frequency) the symmetric so- 
lution becomes unstable at the bifurcation point, while the new asymmetric 
state is nonlinearly stable. We expect that anti-symmetric states (nonlinear ex- 
cited states), corresponding to points on the disconnected branch to the left of 
the symmetric branch, are unstable to generic perturbations due to a nonlinear 
resonance phenomena, and that for large t, the system's energy resides in the 
ground state and radiation modes; see, for example, | 14II15 |. 

2 Overview of Gross-Pitaevskii / NLS 

In appendix^we review the well-posedness theory of 311. 411 . The initial value 
problem is well-posed in C° {R^ ; {R^)) . 

In this section we focus on nonlinear bound states of (|1.3|1 . These are solu- 
tions of the form: 

Mx,t) = e'^'V(a;;A), VGi'- (2.1) 

For the linear Schrodinger equation, the functions are eigenstatcs of a Schrodinger 
operator: —d^-\- v(x) and satisfy the equation 

-V-xx + v{x)ij = -AV- (2.2) 

In the case where v{x) is given by (|1.4|) , with L sufficiently large, the set of bound 
states is two-dimensional, spanned by symmetric and antisymmetric functions. 
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Figure 1: A bifurcation diagram is given for the case e — 1, L = 2. The curve 
on the left bifurcating from M — represents the antisymmetric pulse, whereas 
the curve on the right represents the symmetric pulse. At the threshold value 
A/thresh, an asymmetric pulse bifurcates from the symmetric pulse. At Af— levels 
where the asymmetric pulse is present, the symmetric pulse is unstable. 



For nonlinear bound states satisfy the equation 

+ {vix) - 2\^\'')^ = -A^^. (2.3) 

For the translation-invariant case, w = 0, there is a well-known family of soli- 
tary ("soliton") traveling wave solutions. These are Galilean boosts (see (jA.2|) ') 
of the basic solitary ("soliton") standing wave with hyperbolic secant spatial 
profile. For a non-trivial potential v{x) as in (|1.4|l . the equation is no longer 
translation-invariant and we have "defect" or "pinned" states. These pinned 
nonlinear bound states are plotted in the bifurcation diagram of figure ^ In 
the translation invariant case, v = 0, the nonlinear bound states reduce to the 
branch of NLS solitons, bifurcating from the zero state at zero frequency, the 
endpoint of the continuous spectrum of the linearized operator —d^ about the 
zero state. For v{x) nontrivial, there is a family of defect states, bifurcating 
from the zero state at each linear eigenfrequency of the linearized operator, 
—d^ + v{x) and in the direction of the corresponding eigenfunction; see for 
a general discussion. 

Stability of nonlinear bound states: An important characterization of the non- 
linear bound states is variational. The advantage of the variational characteri- 
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zation is that it can be used to establish nonhnear stabihty of the ground state; 

see [niinj. 

Theorem 1. (I) Nonlinear bound states can be constructed variationally by 
minimizing the Hamiltonian, Ti, subject to fixed L^-norm, J\f: 

min n[<p], Af[ip] ^p. (2.4) 

Such minimizers are called nonlinear ground states, '0Dcf(p)- The associated 
frequency, —X'^{p), arises as a Lagrange multiplier for the constrained varia- 
tional problem (|2.4(l . As p ^ 0, —X^{p) tends to the ground state eigenvalue of 
-dl + v{x). 

(II) Ground states are nonlinearly orbitally Lyapunov stable, i.e., if the ini- 
tial data are close to a soliton (modulo the phase and translation symmetries 
of GP / NLS), then the solution remains close to a soliton in this sense for all 
t S (— oo, oo). 

Although, the above theorem establishes the ground state as a dynamically 
stable state, it leaves unaddressed the question of detailed spatial characteristics 
of the ground state. As we shall now see, dynamical systems methods can, in 
some cases, be used to study this question globally and in detail. 



3 Nonlinear bound states - dynamical systems 
approach 

Dynamical systems methods are particularly well-suited to studying the exis- 
tence and nature of nonlinear bound states for the one-dimensional problem 
p. 311.4(1 . since the equation is piecewise autonomous. In particular, we study 

-i'xx + e v{x)il) - 2ip^ ^ -AV (3.1) 

where v{x) ~ —S{x + L) — 6{x — L). With a quick apology for the following 
non-physical substitutions, we rescale the variables as follows: 

X = z/\ 

V'(x) = A (j){\x) = A (j){z). (3.2) 



In these variables, we have 



,-2 4>''--^v{zl\)(t>. 



A2 

This is equivalent to the following system with matching conditions: 



z 



z ^ ±XL (/)" = - 2 03 



±XL 



0(z+) =0(z-) 

0'(z+) = 0'(z-)^f 0(z), 
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where ' represents differentiation with respect to z and f{z^) = hnip^o f{z±p). 
As usual, we can convert a second order differential equation into a pair of first 
order equations. In particular, if we let (m, v) = (</>, <^'), then we have 



A search for bound states (standing waves) in this system can then be under- 
taken as follows: 

1. Since any standing wave must decay as z ^ — oo, a potential standing 
wave solution must originate (in the phase plane) along the global unstable 
manifold of (0, 0). Without loss of generality, we consider only the portion 
of this manifold that lies in the first quadrant. 

2. When this solution approaches the first well a,t z — —XL^, we may shift z 
so that the point lim^^.;^^- {u(z), v{z)) can be chosen anywhere along the 
unstable manifold. (The value of u{—XL) can be thought of as a shooting 
parameter.) At this chosen point, the solution jumps vertically according 
to the matching conditions for v. 

3. After completing the vertical jump at z = —XL, we again flow the solution 
forward for "time" 2XL according to the ordinary differential equation, 
until we approach the second well at z — +XL. 

4. When z — +XL, the solution again jumps vertically in the phase plane 
according to the matching condition. The phase plane solution that we 
have been following represents a standing wave if and only if this second 
jump lands our solution exactly upon the stable manifold, in accordance 
with the condition that u{z) ^ as z — > +oo. 

For fixed values of the parameters e, A and L, it may happen that there 
arc multiple standing waves present, corresponding to different values of the 
function at the first defect. In accordance with the search method described 
above, such solutions are shown in the phase plane in figure |21 

3.1 The linear limit; infinitesimal solutions 

To illuminate this search for standing waves, it is instructive to first consider 
the case of small amplitude solutions. For infinitesimal solutions, we neglect the 
nonlinear term and we are left with the following linear system 



z ^ ±XL 




z 



±XL 




z ^ ±XL 




z = ±XL 
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Figure 2: Standing waves of different forms are found for e = 1, L = 2 and A = 
0.65. Here we see a symmetric wave, an asymmetric wave and an antisymmetric 
wave. 
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For the linear system, a bounded solution again grows from (0, 0) along 
the unstable subspace, jumps to some transient at the first well, evolves further 
according to the linear system and then must jump exactly to the stable subspace 
at the second well in order to decay to (0, 0) as z ^ oo. For a fixed value of L 
and e in this system, we can determine the values of A for which standing waves 
exist. This is pursued below. 

First, the unstable subspace at (0, 0) is given as the eigenvector of the matrix 




associated with the eigenvalue of positive real part. This is the vector (1,1) and 
so any point in the unstable subspace can be written as (</>, <^). Similarly, the 
stable subspace is spanned by the vector (1,-1). The jump at the first well 
is then J : (0, (p) {(f), (1 — j)(f)- Note that this jump will take us to one of 
the interior transient curves (inside the homoclinic loop in phase space) when 
A > e/2 and to one of the exterior transient curves (outside the homoclinic 
loop) when A < e/2. 




0.1 0.2 0.3 0.4 0.5 0.6 0.7 0.8 0.9 1 



Figure 3: A vs. L for e = 1. Symmetric states exist for any A > e/2, while for 
A < e/2, antisymmetric states exist only for well-separation L > e~^. 

In this infinitesimal amplitude limit, we can find the standing waves of this 
system explicitly. Writing the landing point of the first jump in terms of the 
eigenbasis yields 

(^,(l-l)0)=^(l-^) (1,1) + 0^(1,-1). 
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Because the system is linear in this limit, we ignore the ^-terms from here on. 
Evolving this point forward for time 2XL brings us to the point 



t2AL M 1, 1 - T ) ) = ^2AL 



A 



(M) + ^(i,-i: 



^e2^^ ^ -i^ [e2^^ - e^2AL] g2AL _ ±_ ^g2AL ^ g^2AL 



L [^2^^ _ P-2AL1 2AL _ J_ 

2\\- e j,e 

(3.3) 



Now, the solution we are tracing will represent a standing wave if and only if 
this point is taken to the stable subspace of zero at the second jump. By the 
jump conditions 'At z — ±XL and the fact that points on the unstable manifold 
are of the form (0, —(j)), such points have the form 

J-\W^)=[^,[j-l)cp). (3.4) 

Checking to see if the point (|3.3|) is in this form yields the following equation 

(e - A) (e^^^ - ^ (e^^^ - e'^^^)) - A (e^^^ - ^ (e^^^ + e'^^^)) - 0. 

Some algebraic manipulation shows that this equation is equivalent to 



2A \2X-e 

It may be more natural to consider A as a function of L rather than vice versa; 
this would allow us to determine for which values of A (frequencies — A^) sta- 
tionary waves exist. This can be given implicitly as 

A=|(l±6e-2^^) 

In particular, for A^ to be real and positive, we need the argument of the log to 
be positive and therefore we have the following two solutions 

0<A<./2 ian« = ^ln(^) (3.5) 

e/2 < A < e L,y,n = ^ In (^1^) • (3-6) 

The antisymmetric solutions exist only for L > e^^ because the function ianti(A) 
is monotone increasing on [0, e/2), with minimum value e~^ attained at A = 0. 
The two branches of solutions (symmetric and antisymmetric) are shown in 
figure 01 The results are summarized in the following proposition: 

Proposition 2. In the small amplitude limit, with the wells separated by any 
positive distance L > 1/e, there exists a symmetric standing wave with e/2 < 
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A(L) < e and an antisymmetric wave with < X{L) < e/2. These states are 
unique (up to linear scaling) and no other standing waves exist. Finally, for eL 
sufficiently large, 

A±^|(l±e-^^), A+-A_^ ee-^^ 

3.2 Nonlinear regime; solutions with larger amplitude. 

For solutions with nontrivial amplitude, things get much more interesting (as 
we have seen in our initial numerical search) . Because of the periodicity of the 
transient states, there are many possible solutions for this system (especially for 
larger values of L). We show particularly that in addition to the symmetric pulse 
(and the antisymmetric pulse, for L values at which it exists) described above, 
there is also an asymmetric pulse that bifurcates away from the symmetric pulse 
at a particular value of the boson number A/". (These results follow for a focusing 
nonlinearity, g < 0. In the defocusing case, g < 0, an asymmetric pulse can be 
seen to bifurcate from the antisymmetric solution.) 

For fixed values of the parameters L and e, we will consider what sorts of 
positive-valued standing waves exist for varying values of A. These solutions 
can be determined by adapting the previous search technique in the following 
way: 

1. For this system, the branch of the unstable manifold in the first quadrant 
is given by the explicit expression 

VF" = |(ii,w) : u = wv/T^^, < u < l| . (3.7) 

We now search for standing waves by evolving this entire curve via the 
matching conditions and the evolution equation. 

2. For a particular value of A, application of the matching conditions at the 
first defect yields the curve 

JiWn^{iu,v):v = u(^VT^^j)}; (3.8) 

see the dashed curve in figure El 

3. We then evolve this curve for time 2XL, yielding the curve 

T2XL {J{W^)) ■ 

4. Finally, we compare the evolved curve T2\l ('^(T^")) with the set of points 
that jump to the stable manifold at the second defect 

J-\W') ^{{u,v):v = -u (Vl -u^-j)}. (3.9) 
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This agrees with (13.411 in the linear hmit. Notice that the set J~^(W^) 
is simply the reflection of the set J{W'^) through the v-axis. Each non- 
trivial point of intersection of the two sets T2XL (J(tV)) and J^^{W'') 
represents a standing wave of this system. 



Remark 3. Because the curve J{W") intersects many of the periodic transients 
more than once, it seems plausible that there may be non-symmetric solutions 
for some values of L. We will make that clear by pursuing the strategy above. 

Proposition 4. Let e and L remain fixed. For all A > Agym (where Xsym is the 
value of X for which a symmetric solution exists in the linear limit, as in 1^3.6]) ). 
there is at least one symmetric solution standing wave. 

Proof. This is a simple argument regarding the continuity of the curves J~^{W'^) 
and T2AL {J{W'^)). We consider only the portions of these curves interior to the 
homoclinic orbit of the autonomous system. The portion of J(PF") inside this 
loop has one endpoint at (0, 0) and the other endpoint on the curve W in the 
lower half-plane; see figure 0] 

For A > Asym, the initial portion of the curve T2\l leaving (0,0) is 

between J~^{W'^) and - in particular, above the curve J~^{W^). However, 
the other end must remain on the invariant manifold , and this manifold lies 
below the curve J~^{W^). Therefore, since T2\l (J(VF")) is a continuous curve, 
it must intersect the curve J~^(W^*) in at least one point. Therefore there is at 
least one standing wave in this system for all A > Agym- 

Finally, notice that if there is only one point of intersection, then the standing 
wave corresponding to this intersection point will be a symmetric function of 
position. This follows from the reversibility of the autonomous portion of the 
evolution equation. □ 

Knowing that one or more standing waves exist, we now look into the forms 
that such waves take. In the following, we will look for symmetric standing waves 
that are concentrated at the defects; in particular, we hope to find profiles that 
achieve their maximum values at the two defects. We take a slight change of 
view in the following, considering e and A as fixed and asking for which values 
of the spacing parameter L do such standing waves exist? 

Proposition 5. Let e and A be fixed parameters. There is a symmetric solution 
concentrated at the defects in each of the following regimes: 

1. For I < A < e, there is such a solution for all L > Lsy,n(A) = ^ In ^ 2\-e ) " 

2. For e < A < \/2e, there is such a solution for all values of L. 

3. For A — \/2e, there is such a solution for all L > (tan^^ 2) . 
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Tj, ^ |_(J W") and J"' W= ^ JJ W") and J"' W= 




Figure 4: For all values of A > Asym there is an intersection of T2AL {J{W^)) and 
J~^{W^) - and such an intersection represents a standing wave. As A increases, 
more intersections may arise. 

Proof. We will prove this proposition by considering the amount of "time" it 
takes each point on the curve J{W") (from the point where the curve enters the 
lower half plane until it meets the manifold W^) to reach its reflection on the 
curve J~^{W''). For fixed value 2L, we show there is a point on J{W^) which is 
mapped to its reflection on J~^(M^*) in time 2L. We can define the function L{u) 
as the "time" it takes the solution beginning at a point (m, u{y/l — u-^ — e/^)) 
(meeting the previous parenthetical requirement) to reach the symmetric point 
{u,-u{VT^^-e/X)). 

1. The result follows quickly in the regime | < A < e, by noting that as 
M 0+, L{u) — > Ls as given in the earlier proposition. Also, as w ^ 

(^^yT~~(eJ2X)^^ , L{u) 00. A continuity argument completes this 
case. 

2. For A > e the curve J{W") leaves the origin into the upper half plane. It 
reaches the u-axis when u = -^1 — (e/A)^ and reaches the stable manifold 

when u = y^l — (e/2A)^. In this regime, it is clear that L{u) ^ as 

u (^1 - (e/A)2^ ^ and L{u) ^ cx) as w ^ (^y/l - (e/2A)2^ 

3. When A = V^e, the curve J(VF") crosses the u-axis at the elliptic fixed 
point (l/-\/2, 0). A careful study of the linearization near this fixed point 
reveals the lower limit given above. 

□ 

Now we note an important consideration concerning the form of these sym- 
metric pulses. In particular, for fixed values of the parameters e and A, there 
is a distinguished value of the parameter L so that the curve J{W") is tangent 
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to the symmetric standing wave as we approach the first defect from the right. 
The value of the standing wave at each defect can be quickly computed for this 
threshold value, 



12-(i)^-lVl2+(f)^ 
18 

and the threshold value of L, Lthrcshi can then be computed simply by noting the 
"time" it takes to flow from (uthrosh, Wthiosh) to its reflection (wthrosh, -^thresh)- 
In figureEl we show the energy of the symmetric standing wave at the thresh- 
old values of L. We will now show that beyond this threshold, there is a pair 
of asymmetric standing waves (in addition to the symmetric and antisymmetric 
waves discussed previously). 



'^thresh 




Stable symmetric 
ground states 



Figure 5: Threshold values of the particle number Athrcsh for varying values of 
L are shown. At these threshold values, there is a bifurcation of asymmetric 
pulses (see Proposition inj and beyond these values, the symmetric standing 
wave will be unstable (see Theorem (HJ. The point along the curve marked with 
a solid circle corresponds to the special value A = V^e, where the point Mthrosh 
is the flxed point in the phase plane. Along the curve to the right of this point, 
the symmetric standing waves at threshold have the desired two peak profile; to 
the left, the symmetric standing waves at threshold have a single peak centered 
between the two wells. 



Proposition 6. Let e and A remain fixed. For L > Lthicsh(A), there exists a 
pair of non-symmetric standing waves. These waves leave the symmetric wave 
in a pitchfork bifurcation at L = Lthresh- 
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The solution u(z) „ Symmetric and asymmetric waves u(z) 

0.26 1 1 1 1 1 1 1 1 1 1 1 ■ I ' ' ' ' ^ 




Figure 6: The solutions corresponding to the various intersection points from 
the phase plane diagram in figure 0] are shown here. 

Proof. Consider e and A fixed and positive, with L so that L > Lthrcsh(A). 
For this value of L, the symmetric standing wave has a value (Mo,fo), with 
uo > Mthrcsh, approaching the first defect from the right. This wave then must 
pass through J~^{W'') in the phase plane an additional time before the reaching 
the terminal point at the second defect. 

In the phase space, the part of the standing wave between the two defects 
can be written as the curve 

To = {T,{{uo,vo)):0<z<2XL}. 

Note that the curve Tq simply traces out the portion of the periodic orbit con- 
necting the symmetric points (mo,wo) and {uq, —vq). 

We consider an arbitrary parameterization (ut,Vt), t g [0, 1] of the portion 
of the curve J(VF") between {uo,vo) and its nontrivial intersection with the 
stable manifold W^, so that (uo,wo) is as above and {ui,vi) € W. Then we 
can define 

= {T, {{ur,vr)) : < z < 2XL} , 

for all r G [0,1]. Notice that the curve J~^{W^) is transverse to all of the 
periodic orbits in the phase plane away from the fixed point (l/-\/2, 0). Therefore 
if a curve Tr crosses J~^{W), it must do so transversely. 

Recall that the curve Fq crosses J~^{W'^) prior to its terminal point at 
the second defect. (This earlier crossing does not represent a standing wave, 
of course, since it occurs for a value of z not corresponding to a defect.) On 
the other hand, for r sufficiently close to 1, the curve F^ no longer crosses 
J~^{W^). Continuity considerations therefore dictate that there must be a 
value of r between and 1 so that T2XL {{ut,Vt)) hes on the curve J~^(PF'*). 
This value clearly does not represent a symmetric state since uq > Wthrosh, while 
the first point of intersection is smaller than Mthresh- 

There is a second non-symmetric standing wave that, in the phase plane, 
has the form of the wave just found refiected across the u-axis. 
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This pair of non-symmetric solutions is originally produced in a pitchfork 
bifurcation when L = Lthrosh ~ at this parameter value the transversality of 
J~^(W*) and the periodic orbit carrying the symmetric standing wave is lost. 

□ 



4 The eigenvalue problem and stability. 

In the previous section, we found various standing wave solutions of the equation 

ii/'t = -i^xx - 2|V'PV' + e v[x)ip 

where v{x) is the double well potential —6{x — L) — 6{x + L). In this section, 
we will analyze the stability of the various branches of solutions constructed, 
by studying the eigenvalue equation, associated with linearization about a fixed 
nonlinear bound state. 

First we change coordinates into the frame of the standing wave 

ipix.t) = e'^'V(a;,i) 
and arrive at the equivalent equation 

^<t>t = -4>xx + [>? - (2|<^P + e v{x))] 4>. (4.1) 
We write this equation in real and imaginary parts (<j) — v + iw) 



wt = Vxx + [2{v'^ + w'^) + e v{x)]v - 



(4.2) 



and formally consider small perturbations of the standing wave. Since the stand- 
ing wave (let's call it u{x)) is real-valued, we then write (v{x,t),w{x,t)) = 
{u{x),Q) + 5 {p{x,t),q{x,t)) with 5 small. Leading order in 5, yields the lin- 
earized dynamics 



Setting 



Pt = -qxx - [2u'^ + e v{x)]q + A^g, 

qt = Pxx + [2u'^ + e v{x)]p - iu'^p - A^p. 



+ [2u^ + e v(x)] - 



L+ = A_ + [2u^ + ev{x)]+4u^-X^, 
equation (|4.2() can be rewritten as the system 

P \ _ f -L_ \ f P \ _ f P 



q Jt \ L+ J \ q J \ q 

Of interest is the spectrmn of this operator TV. For the ground state (see section 
El), the spectrum of N lies entirely on the imaginary axis and modulo zero modes. 
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associated with symmetries, exp(tA^) is bounded on and indeed, the solution 
is nonUnearly stable |17[ [T5| . 

A state is unstable if one can exhibit an eigenvalue of the operator N with 
positive real part. Nonlinear instability is a consequence; see the appendix of 
Unfortunately, even though the spectrum of L+ and L_ can, in principle, 
be determined using Sturm-Liouville theory, it is not obvious how to construct 
the spectrum of N using this information. However, this is exactly the operator 
studied by Jones in in an application to optical pulses in nonlinear waveg- 
uides. In that work, Jones interprets the search for a real positive eigenvalue 
as a shooting problem in the space of Lagrangian planes. Using a winding ar- 
gument in this space, he derives a relationship between the number of positive 
eigenvalues of the simpler operators and L_ and the number of positive 
eigenvalues of N. Let 

P = number of positive eigenvalues of and 
Q = number of positive eigenvalues of 

Then we have 0111: 

Theorem 7. If P — Q ^ 0, 1, there is a real positive eigenvalue of the operator 
N. 

An alternative, PDE approach to this result was developed by Grillakis [Sj. 

From Sturm-Liouville theory, P and Q can be determined by considering the 
behavior of solutions of L^p = and L^q = 0, respectively. In particular, P 
and Q equal the number of nodes of the associated solutions p and q. Since 
L^q = is satisfied by the standing wave u{x) itself, it quickly follows that 

Q — number of zeros of the standing wave u{x). 

For all positive-valued pulses, then, Q = 0. 

So we have that if P > 1, then N has a positive eigenvalue and the standing 
wave u{x) is unstable. The operator L_|_ acts as the equation of variations for 
H4.1|l . and a natural interpretation of the equation of variations is that it carries 
tangent vectors under the flow. So if we initialize q{x) as a vector tangent to 
the unstable manifold at the point u{x), then we can calculate P by counting 
the nodes of the resulting solution q{x). 

In the context of vectors, a node of q{x) is a vector in the phase plane 
pointing straight up or down. And so, P can be computed by counting the 
number of times that the vector initialized tangent to the unstable manifold 
passes through verticality. In particular, considering only the angle of this vector 
and not its magnitude, P corresponds to the largest number of integer multiples 
of 77 through which this vector rotates as it is evolved under the linearized flow. 

The only difficulty in the computation is this: because of the jumps at 
the defects, i.e., the inhomogeneity of the evolution equation, a vector q{x) 
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initialized tangent to the nonlinear bound state solution will not remain tangent 
to this solution once we reach the first defect. 

However, using the simple matching conditions for these vectors at the wells, 
we have the following theorem: 

Theorem 8. For L < ithrcshj the symmetric solution has P < 1. For L > 

ithrcshj the symmetric solution has P > 2. Therefore, after the bifurcation of 
the asymmetric waves, the symmetric wave is unstable. 

Proof. The result follows simply by evolving a vector (initialized as a tangent 
vector to the upper branch of the homoclinic) around each standing wave under 
the influence of the equation of variations. This can be visualized by noting that, 
at the first defect, the tangent vector to the unstable manifold is taken to 
a vector that is tangent to the curve J(VF"). The observation of instability is 
then reduced to simply noting the manner in which this landing curve J(W") 
intersects the periodic curve in the phase plane representing the standing wave 
between the defects. 

1. At the threshold value L — ithrosh itself, the curve J(VF") is tangent to 
the transient to which the solution jumps at the first well. The varia- 
tional solution, then, will still be tangent to the transient curve (although 
pointing backwards). Hence, we can continue to trace the solution exactly 
by following the tangent vector around this transient to the second jump. 
Then, as symmetry dictates, this vector is taken to a tangent vector to 

at the second jump. This vector is seen to pass through verticality 
only once (at t = 0) and hence has P = 1. 

2. For L < Lthrcsh, the vector jumps through a smaller angle at each well 
than it does at the bifurcation point itself. At the first defect, it then 
points along the curve J(W^"), which in this case points up and into the 
internal ellipse. Because solutions are unique, the tangent vector to this 
ellipse provides a bound for our solution as it evolves, and the vector in 
question must retain its orientation with respect to the tangent to the 
ellipse; in particular, it rotates clockwise and must be pointing "below" 
the ellipse at the second defect. However, the vector that is taken to the 
curve is the vector tangent to J~^(VF^), and comparing these two 
vectors we quickly see that P < 1 in this case. 

3. For L < Lthresh, the vector jumps through a larger angle at the first defect 
and now points along J{W") down and out of the internal ellipse, and so 
the tangents to the ellipse now provide a bound on the other side. At the 
second well, the solution must still point out of the ellipse and comparison 
with the vector tangent to J^^{W^) forces P > 2 in this case. 

These three cases are sketched in figured □ 
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Figure 7: The evolution of the tangent vector is shown for solutions with varying 
values of L. In the first figure L — ithrcsh and the evolved vector is parallel to 
the solution throughout its evolution. In the second figure L < Lthrcsh and the 
vector moves through a smaller angle at each defect and we see P < 1. When 
L > ithrcsh as in the third figure, the vector moves through a larger angle at 
each defect and forces P > 2. 
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A Wellposedness of Gross-Pitaevskii / NLS 



Structural properties of NLS: The Gross-Pitaevskii / NLS equation is a Hamil- 
tonian system, whieh can be written in the form: 

iut = —Hgp[u,u], (A.l) 
ou 

where Ti.Qp[u, u] denotes the Hamiltonian (|1.5(l . Invariance with respect to time- 
translations imphes that Hgp[u, u] is conserved by the flow generated by ljA.l|l . 
Additionally, invariance under the transformation u e^^u, ^ G i?^ implies 
that A/'[m, u] , the norm defined in , is conserved. 

For the spatially translation invariant case, NLS has the Galilean invariance: 

u{x,t) ^ u(a;-si,t) e'^^"-^"'*), s e R\ (A.2) 

Well-posedness theory: We present a sketch of the existence and uniqueness of 
solutions to (|1.3(l : see also The functional Hgp[-] and N[-] are well defined 
on H^{R^), the space of functions /, for which / and dxf are square integrable. 
It is therefore natural to construct the flow for initial data of class H^. In fact, 
it can be shown that, for initial conditions uq — u{x,t = 0) G H^{R^), there 
exists a unique global solution of NLS, u G C^{R^\ H^{R)), in the sense of the 
equivalent integral equation: 

u{t) = U{t)ua +2i [ U{t- s) \u{s)\'^u{s) ds, (A.3a) 

U{t) = exp(-ifft), H = -dl + v{x). (A.3b) 

The spectral decomposition of H is known explicitly and can be used to 
construct U{t) expHcitly. 

To show the existence of a solution to (|A.3a|l in H^, we must show the 
existence of a C'^ {R^ ; (R)) fixed point of the mapping u{x,t) J[u]{x,t), 
given by the right hand side of (|A.3a() . 

We now outline the key ingredients of the proof. To bound J[u] and its first 
derivative in L^, we introduce the operator A = I + PcH, where Pc denotes the 
projection onto the continuous spectral part of H. Note that A is a. nonnegative 
operator, since the continuous spectrum of H is the nonnegative real half-line. 
Moreover, we expect ||-42/|| ~ ||/||_ffi = ||(^^ d^)'^f\\L^- In fact, we shall use 
that the following operators are bounded from to L^: 

AHl-d'J-i, A~i{I-d'J. 

This follows from the boundedness of the wave operators on |16|. Therefore, 
we have an equivalence of norms 

Ci WfWm < WAifh^ < C2 WfWm. (A.4) 

Our formulation (|A.3ap and introduction of A are related to the nice property 
that A^ and hence also functions of A, commute with the propagator exp(— i_ff<). 
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We shall also use the Sobolev inequality: 

\f{x)\' < C\\f\\L2\\d,f\\L2, (A.5) 

and the Leibniz rule |1U| : 

(A.6) 

Since U{t) is unitary in L^, we have 

\\Ah[u]\\L2 

< M^uoIIl^ +7 / \\A^'\u{s)\^u{s)\\l2 ds 

JQ 

= WA^mw^. + jl\\[AHi-di:r^)-{{i-d',)i\u{s)\M^))\\L2ds. 

(A.7) 

By (TOll . (EHJ, and ifTfill . 

PNWIIhi < C\\AJ[u]\\l2 < CiWuaWm + C2T sup \\u{s)\\%,. (A.8) 

se[o,T] 

Now assume that u is such that sup^gjo ||w(s)||//i < 2Ci. Then, by HA.8|) . 
choosing T <Ti sufficiently small, supsg[Q j^j || J[u](s)||/fi < 2Ci||uo||//i and so 

\\J[um\\m < C.WuoWm + C2T {2Ci\\uo\\mf . 

It follows that for < T < Ti sufficiently small, the transformation J[-] maps a 
ball C'^{[0,T]; H^{R)) into itself. A similar calculation shows that 

\\J[u]it) - J[v]mm < KT {2Ci\\u4H^f sup \\u[s) - v{s)\\m , (A.9) 

se[o,T] 

and therefore for < T < T2 < Ti, the transformation J[-] is a contraction 
on this bah. Therefore, J[] has a unique fixed point in C°([0, T]; for T 

sufficiently small and local existence in time of the fiow follows. Global existence 
in time follows from the d priori bound on the norm of the solution implied 
by the time-invariance of norm. A/", and the Hamiltonian, TYgp- 
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